Non-universal dynamics of staggered non-equilibrium 
particle systems and Ising chains 

R. B. Stinchcombe 1 , J. E. Santos 1 and M. D. Grynberg 2 

1- Department of Physics - University of Oxford 
Theoretical Physics, 1 Keble Road, Oxford 0X1 3NP, UK 
2-Departamento de Fisica, Universidad Nacional de La Plata. 
(1900) La Plata, Argentina 
PACS numbers: 05.70.Ln, 75.10.Hk, 02.50.Ey 



Published in J. Phys. A 31, 541-549 (1998). 



Abstract 

Non-universal dynamics is shown to occur in a one-dimensional non-equilibrium system of hard-core 
particles. The stochastic processes included are pair creation and annihilation (with rates e and e') and 
symmetric hopping rates which alternate from one bond to the next (pa, Pb)- A dynamical scaling 
relation between the relaxation time and the correlation length in the steady state is derived in a simple 
way for the case e' > pa > Pb » e. We find that the dynamical exponent takes the non-universal value 
* = 21n(£)/ln(fl£). 

For the special condition e + e = pa + Pb, where the stochastic system is in principle soluble by 
reduction to a free fermion system, the model is mapped to the Glauber dynamics of an Ising chain 
with alternating ferromagnetic bonds of values Ji and J2, in contact with a quantum thermal bath. The 
full time dependence of the space-dependent magnetisation and of the equal time spin-spin correlation 
function are studied by writing the master equation for this system in the quantum Hamiltonian formal- 
ism. In particular we obtain the dispersion relations and rigorously confirm the results obtained for the 
correlation length and for the dynamical exponent. 
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1 Introduction 

This paper discusses non-universal critical dynamics of non-equilibrium particle systems, in which dynamic 
exponents depend on microscopic parameters (ratios of transition rates). 

Non-equilibrium particle systems with stochastic dynamics show properties ranging from steady state 
phase transitions to turbulence, shocks, and non-equilibrium analogues of strong fluctuations and critical be- 
haviour jj], H . As with equilibrium critical phenomena, associated critical exponents are typically dependent 
on symmetries and dimensionalities and (for dynamics) conservation laws, but not on microscopic details. In 
particular the dynamic exponent z in many simple systems describes diffusive relaxation towards equilibrium 
(z = 2), while in others it takes other parameter independent values, e.g. z — 3/2 in the one-dimensional 
noisy KPZ or Burgers equations ||, [|, ||, 

Simple hard-core particle models of non-equilibrium stochastic dynamics can be mapped to quantum sys- 
tems and in particular to spin models. This latter mapping is achieved either by a pseudo-spin representation 
of the particles or by associating particles with domain walls |[ [| |l0[ [H] ■ 

But it is well known that non-universal behaviour can exist in certain non-uniform spin systems with 
Glauber dynamics jl^, such as that providing a simple model for freezing into non-equilibrium states and 
glassy dynamics [[l3[ [l5| Here the relationship r ~ £ z between the relaxation time r and the 

equilibrium correlation length involves a parameter-dependent dynamic exponent z |l7[ [l8[ [l9f . This suggests 
that analogous non-universal behaviour can exist in non-uniform stochastic particle systems, and this is 
demonstrated here. 

The model considered is defined on a one-dimensional lattice with L sites. We divide this into two sub- 
lattices, A and B, containing even and odd sites respectively. Each site can be occupied by at most one 
particle. A pair of particles in adjacent sites can annihilate with rate e'. A pair can be created on two empty 
adjacent sites with rate e. A particle can hop to an empty nearest neighbour site with rate pa if it is in 
sub-lattice A and with rate pb if it is in sub-lattice B. 

The correlation lengths characterising the steady state particle separations on each sub-lattice can be 
obtained by dynamic balance conditions within a mean field approximation. The characteristic time for 
diffusion across the smaller of the lengths, £, can then be estimated from the hopping process. This is 
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particularly simple to do when e.g. the annihilation rate is much greater than the creation rate (when £ 
is large) and one hopping rate is much greater than the other. This simple procedure already provides the 
non- universal behaviour of the dynamic exponent z. 

For a special case where the rates satisfy one constraint, the quantum spin Hamiltonian representing the 
particle dynamics reduces to a free fermion form after a Jordan- Wigner transformation |20|] . The alternative 
and simpler way used here of achieving an exact solution for this case is by mapping the stochastic particle 
system to the Glauber dynamics of an Ising model with alternating bonds J\ , J2 jlO| ^l) . 

The correlation lengths and relaxation time readily provided by the simple argument in the particle 
picture thus translate, for the special rate relation, into results for the alternating bond Ising model, where 
they can be confirmed by exact calculation on the spin model. 

The Glauber dynamics of this model is that resulting from coupling the static Ising system to a quantum 
thermal bath and considering the Van Hove limit of a weak interaction and very large times. By writing 
the master equation which describes the Glauber dynamics in a quantum Hamiltonian formalism, we can 
easily obtain the equation determining the time evolution of the equal-time spin-spin correlation function. 
From this quantity one can compute the equilibrium domain wall density in the Glauber problem which is 
equivalent to the steady state density of particles in the stochastic particle dynamics Jl(| . By taking 
the Fourier transform of this equation we obtain a system of four linear equations which determine the 
dependence of the eigen-frequencies with the wave vectors of the two spin excitations. The solution of the 
associated secular equation gives us the dispersion relation. The behaviour of the frequency at low wave 
vectors determines the critical exponent z. We fully confirm the results obtained previously. 

The structure of this paper is as follows: in section 2, we define the dynamics of the particle system 
in terms of the constituent stochastic processes (diffusion, annihilation and creation). Using the simple 
mean-field and random-walk arguments referred to above we obtain the relation between the relaxation time 
and the equilibrium correlation length and hence we determine the critical exponent z. In section 3 we 
present the Glauber dynamics of the alternating-bond Ising system, giving the rates of transition between 
the different configurations. We then go on to relate the spin flip processes to domain wall hopping, creation 
and annihilation, which translate to the particle processes. In section 4, we write down the master equation 
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in a quantum Hamiltonian formalism and derive the equation giving the time evolution of the spin-spin 
correlation function. Hence, we obtain the dispersion relation and from it we can again compute z. Finally, 
in section 5 we present our conclusions. 

2 The reaction-diffusion system and its non-universal dynamics 

The hard-core particle dynamics is as specified in the Introduction: a particle on the even sub-lattice may 
hop to its left nearest neighbour site, provided it is empty, or to its right nearest neighbour (if empty) at 
rate pa] the corresponding rate for hopping from the odd sub-lattice is pb] in addition, pairs of particles can 
appear at (or annihilate from) adjacent empty (or full) sites at rate e (or e'). These processes are depicted 
in Table 1 (see below). 

We now consider the description of the steady state. Here the particle densities pa and ps on the two 
sub-lattices will be uniform. Since the system is in equilibrium the processes of creation and annihilation 
have to balance each other. So one concludes that provided the mean field approximation applies, e' paPb = 
e (1 — pa){1 — Pb)- The same must be true for the processes of diffusion. Hence pa pa(^ — Pb) = Pb (1 — Pa)Pb- 
One obtains from these relations p* ~ , X , A , where x& — f-^ 2 -) 1 ^ 2 and on ~ = XB , where xn = (A^) 1 ^ 2 . 
These densities allow us to extract two distinct separation lengths, the average distances p~ x between particles 
on each sub-lattice a. 

We will be particularly interested in the critical situation where these two characteristic lengths are 
both large. That occurs when e' 3> e, giving large particle separations, or where e ^> e' (large vacancies 
separations). For convenience we consider only the former case. It is also convenient to take pa 3> Ps, 
to widely separate the two lengths. Both are still large provided e'/e 3> Pa/pb- The shorter, controlling 
correlation length £ is then 

i - Pb 1 - ^B 1 = {z'PB/ePA) 1/2 ■ (1) 

Now to determine the rate of approach to this steady state for e' 3> e one uses a modification, appropriate 
to the particle dynamics, of an argument given by [53]. For large times the evolution of the system is deter- 
mined by the limiting relaxational process, namely particles diffusing until they meet, when they annihilate. 
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For situations close to the steady state, i.e. for long time critical dynamics (starting from generic non- 
equilibrium initial states), the characteristic time involved is that for particles to diffuse across the shorter 
steady state separation length, £, and to annihilate. Now to traverse two adjacent bonds both rates pa and 
Pb enter and the effective diffusion rate is ApAPB , Thus the characteristic time is r ~ (e'pn/epA) PA+PB oc £ z 
giving 

z = 21n(-)/ln(^) (2) 
e tPA 

in the regime of validity (e'/e 3> pa/pb 3> !)■ This is a non-universal result, being dependent on ratios of 
rates. The above reasoning is based on the mean field argument that we have used to determine pa and 
Pb- Mean field theory determines correctly the characteristics of the steady state but it fails to predict the 
correct approach to the steady state [jll). This is due to the diffusion aspects of the problem, which are 
relevant in one dimension. The random- walk argument that we have given to determine r incorporates these 
aspects. We will see in the next sections that the same results are obtained by an exact calculation based 
on an equivalent Glauber-lsing dynamics. That equivalence applies when the rates satisfy 

e + e' = Pa+Pb ■ (3) 

This relation is also sufficient to reduce the quantum Hamiltonian representing the particle dynamics to a 
free fermion form, which gives another possible method for solving the particle dynamics exactly. 

Finally, we note that the exponent z changes continuously with the value of the rates of the diffusion 
and reaction processes, as given by equation (|^). A similar example of a critical exponent which changes 
continuously was found in the study of the dynamics of the g-state Potts model J2^| . In this system, the 
fraction of spins which never flip up to time t, r(g, t), decays like a power law r(q, t) ~ t~ s ^ q ' when the initial 
condition is random. The exponent 9{q) varies continuously with q. 

3 An Ising system with generalised Glauber dynamics, and its 
relationship to reaction-diffusion processes 

The Ising system considered is a spin | chain, where a spin on the even sub-lattice is coupled to its left nearest 
neighbour by a ferromagnetic Ising interaction of strength J\ and coupled to its right nearest neighbour by an 
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interaction of strength J2 (Ji > J2) p7p . One can therefore write the Hamiltonian for this alternating-bond 
system in the form 

L 

H = -J2 Jimm+i (4) 
1=1 

where 07 denotes the eigenvalues of af, Ji — J\ for I odd and J; = J2 for I even. The transitions between 
different configurations of the system, i.e. between different sets of values of the spins, are of the Glauber 
type, i.e. the system evolves by single spin flips [ fL2|| . The probability per unit time that a spin will flip from 
its value to the opposite one is taken to be JI?]] 

r 1 

= 2 t 1 ~ 2 + Ti (5) 

with 

7 ± = tanh(AVi + K t ) ± tanh(AVi - (6) 

where X; = Ji/ksT. Notice that if we take J\ = J2 we will recover the usual Glauber rates. It can be easily 
verified that these rates satisfy detailed balance which is a sufficient condition for the steady state distribution 
of the Ising system to be a Gibbs distribution. This dynamics can be derived pit by coupling the Ising 
system, described by Hamiltonian (|]) to an ensemble of free fermion baths with a grand canonical probability 
distribution, i.e. with a density matrix given (in the case of no interaction) by p = exp[— (3(H — fj,N)]/Z 
where H and N are respectively, the Hamiltonian and the particle number operator of the free fermion bath, 
P is the inverse temperature and fi the chemical potential. The quantity Z is the grand partition function 
and is just a normalisation factor. If the particle number of the fermion system is kept fixed, then the 
chemical potential is, for low enough temperatures, essentially equal to the Fermi energy of the system. The 
coupling is done via an interaction operator which couples the a x component of each spin to a thermal bath 
whose probability distribution is the one given above. The different copies of the thermal bath are totally 
uncorrelated. If one takes the joint limit of the interaction strength A going to zero and t — > 00 such that 
t A 2 is a constant (limit of Van-Hove) J25|, ^(| one obtains the transition rates given by (j|) , in the condition 
that the chemical potential (Fermi energy) of the fermion bath is much larger than the couplings Ji and 
J 2- If that is not the case then the transition rates still have the form (j|), but the parameter T is no longer 
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Initial State 


Final State 


Transition Rate 


Totot 

1 1 1 


T • 1 • T 


E [1 - tanhfKi + K 2 )} = e 


T • 1 • T 


Totot 

1 1 1 


E [1 + tanhfKi + K 2 )l = e' 


total 

1 1 •4' 


t • 1 o 1 


| [1 + tanh(Ki - K 2 )] = Pa 


T« W 


T°W 


I [1 - tanh(K 1 - K 2 )] = p B 



Table 1: Processes of transition for a spin in the even sub-lattice and the equivalent particle processes. The 
spin which is to flip is the central one. The coupling strength of the left bond is J\, and of right bond Ji- 
An empty circle in the dual lattice is to be identified with a vacancy and a full circle with the presence of a 
particle. For a spin in the odd sub-lattice the rates of the last two processes should be interchanged. 



a constant and depends on the value of the neighbouring spins. We will consider that we can take T as a 
constant. 

It is well known that there is a duality relation between Glauber dynamics and the reaction-diffusion 
model of hard-core particles discussed above . To see this we take a given configuration of the Ising 

system and consider the lattice of sites located in the middle of the bonds between the Ising spins (dual 
lattice). If the neighbouring Ising spins have different signs then we place a particle at that site of the 
dual lattice. Otherwise we leave the site empty. In this way we map domain walls in the Ising system to 
particles in the dual lattice. It can be shown that this mapping has a precise mathematical meaning |27| . 
The different possible processes of transition for a given spin, its translation in terms of particle processes, 
the rates associated and their relation are given in table 1. Notice that the sum of the rates of the first two 
processes is equal to the sum of the last two in agreement with (^). This shows that this system is in the 
class of systems that are integrable (the precise meaning of this word will be made clear below) through 
free fermions 0. This allows the computation of correlation functions for a set of distributions of initial 



configurations |11, gTj. Here we will not pursue this point. If we now take the limit of low temperatures 
(K\ ^> K2 1) we see that the rates of pair annihilation (e') and of diffusion from the strong bond to the 
weak bond pa are equal to T. On the other hand the rate of diffusion from a weak bond to a strong bond is 



Pb ~ rexp(— 2(K% —Kv)). Finally the rate of pair creation is e ~ rcxp(— 2(K\ +^2))- So we see that these 
two last processes are exponentially suppressed at low temperatures. Nevertheless they have to be taken 
into account for a proper description of the steady state and diffusive dynamics, as we saw in section 2. We 
will see in the next section that the results given in section 2 may be obtained by an exact calculation. 

4 The master equation in the quantum Hamiltonian formalism 

In this section we will study the time evolution of the spin-spin correlation function. For that we need the 
full master equation for the evolution of the Ising system with rates given by (|5|) . Since the master equation 
is a linear equation a particularly convenient way to write it is to use an operator formalism which assigns 
to each configuration of Ising spins a vector n) in an Hilbert space. The probability distribution for the 
different configurations at a given time t can then be written as a state- vector \^t) = ^2 n P(n,t) \ n) 
[ p8| , where P(n,t) is the probability to find configuration n at time t and is a solution of the master 
equation. The set of different n) is supposed to be orthonormal and complete. One can then write the 
master equation in the compact form dt | W t ) = — T \ } > where T is a linear and in general non-hermitian 
operator that for two state systems with local interactions can be written as a quantum spin Hamiltonian. 
The average values of quantities like the Ising spins can also be conveniently represented in this language as 
(<Jj(t)) = (s I ct| I }, where <rj is the Pauli spin matrix at site j and ( s = J^n I- Substituting the 
state vectors and the operators by their definitions and using the orthonormality relations one sees that one 
obtains the usual definition of average over configurations. One also can, given the correspondence with the 
Schrodinger equation, define an Heisenberg representation of the operators by A(t) = e Tt Ae~ Tt , where A 
is a generic operator. These operators obey the equation of motion 4^ = [T , A(t) ]. It should be stressed 
that this is only a convenient way to represent the master equation for this system and is not related to the 
intrinsic quantum dynamics of the Ising system which is given by the Hamiltonian (Jij). In fact, the dynamics 
is in this case generated by the interaction of the Ising spins with the thermal bath and cannot be deduced 
from the form of (^) alone. For the system that we have studied in the previous section the operator T has 
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the form 

L 

2 



T = ^E^-^^-S^-i+^i)] ( 7 ) 
1=1 

as can be seen if we consider the matrix elements (n' \ T\n) and ( n \ T\ n ) , where n' is a configuration 
differing from n by the flip of a single spin tr/ , say. We obtain 

r l 

(n'\T\n) = -- [1 - -<ji (7;tifTi_i + Ti+i^+i) ] 



r i 

<n|T|n) = - ^[1 - -<Ti (j^cri-i + jf +1 a i+1 )] (8) 
i=i 

where ai, etc., are the values of the Ising spins in configuration n. The first relation gives (up to a minus sign) 
the transition rate between configurations n and n'. The second relation gives the total rate of transition 
out of configuration n. This agrees with the general definition of the f operator p£| ] . The form (Q) and 
associated equations of motion are particularly convenient for they allow the calculation of multiple time spin- 
spin correlation functions. In the same way one can represent the master equation for the particle dynamics 
in terms of a T operator which is also a quantum spin Hamiltonian. The duality relation between these two 
models can be given a precise mathematical meaning by means of a similarity transformation between the 
operators of the two models |27j ). In particular, the T operator of the generalised Glauber dynamics maps 
to the T operator of the particle system. This later operator can be expressed in terms of free fermions by 
a Jordan- Wigner transformation |2p|] and therefore the system is completely integrable. Many correlations 
functions relevant for the study of the particle dynamics and the associate Glauber problem can hence be 
obtained. 

If one now considers the Heisenberg equation of motion for <tJ (t) and takes the average value on some 
arbitrary initial state VP } in the way indicated above one obtains the differential equation giving the time 
development of the average space dependent magnetisation [ fl7| . This equation can be solved by defining the 
Fourier transforms and mj7 of the magnetisation at even and odd sites respectively, with the wave vector 
k equal to k = ^-n with n an integer and such that k lies in the limits — ~ < k < J (we have chosen the 
lattice spacing to be one). Notice that the Brillouin zone has been reduced to half of the size it would have 
if Ji = Ji since the lattice has now a periodicity of two. This procedure gives us a system of two equations 
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coupling and m k for every k. The relaxational eigenvalues of this system are 

ujf = 1 ± [cos 2 k tanh 2 (iTi + K 2 ) + sin 2 k tanh 2 ( J ftT 1 - K 2 )} 1/2 (9) 

where tu^ stands for the lower (acoustical) and upper (optical) branches of this dispersion relation. The late 
(critical) dynamics is determined by the low k modes of the acoustical branch lo^ . The characteristic length 
associated with the decay of the local magnetisation in the steady state can be found from the value of | k | 
(= 2p) making the right hand side of (^) vanish. We therefore take the analytic continuation of @ to the 
complex plane. We obtain the space decaying mode with infinite lifetime (which decays over a characteristic 
length £). For low temperatures £ is given by 

£ ~ exp(2K 2 ) . (10) 

Therefore at low k and low temperatures one can write the dispersion relation as 

oj- = 2rr 2 (i + (fcO 2 ) (ii) 

with £ given by ( pX| ) and 

z = 2 + ±_Jl (12) 

a result given in Jl^| . Now to determine the time evolution of the density of particles we need to compute 
the equivalent quantity in the Glauber language. This quantity is the domain wall local density (Tcfl which 
is given by ^(1 — {ai{t)ai+\{t))). We therefore need the equation giving the time evolution of the spin-spin 
correlation function. Notice that these correlations do not involve only eigenmodes of (0) orthogonal to 
the slowly decaying ones, so their relaxation follows closely the slowest processes in the system. Using the 
transfer matrix formulation we can write ^(s |<r ; z (t)ir^ l (f)| ) = (s \ [T , af (tya^it) ] \ ^ } where | VP ) is as 
above a general initial state. Substituting T by its expression given in ^ we obtain the equations 

d r 

-{cJi{t)cr m {t)) = -2T{ai{t)a m {t)) + -[ 7 +(c7 i _ 1 (t)(7 m (i)} + 7 f {<r l+1 (t)a m (t)) 

+ >y+{<n{t)* m -i(t)) +7mfa(*)°m+x(t)>] (13) 

if I ^ m and the boundary condition (ai(t)ai(t)) = 1. Since these equations are linear this boundary 
condition may be ignored. It can be later enforced by superposing different solutions ||l2|| . To solve these 
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equations we introduce Fourier transforms of the quantities (cri(t)a m (t)) when I and m are even, I is even 
and m is odd, etc. The two wave vectors k and k which now appear are defined as above. We obtain a 
system of four linear equations. The associated secular equation has four solutions corresponding to having 
two excitations in the acoustic branch, one in the acoustic and one in the optical branch (two solutions) and 
two excitations in the optical branch. The frequency Ukk' is given by 

Wfefe' = u k + uj k > (14) 

where uJk and uiy are given by one of the solutions in (|^). If we now take the translational average of the 
correlation function (<7;(t)<7; + i(t)) we see that only the terms for which k' — —k remain in the Fourier sum. 
Hence the decay of density of domain walls (density of particles) is determined by the acoustical branch of 
u>k,—k for l° w k. At low temperatures we just obtain uJk,-k = 2cj,7 with lu^ given ([ll]) . The interparticle 
distance in the steady state is given by the value of k that makes uJk,-k vanish. This coincides with (^p|). 
The relaxation time r is given by r = £ z /4r (t" 1 = uq,o). This translates into (g) using the relationships 
between rates in the Ising and particle pictures (see Table 1) in the low temperature limit. So we recover 
the results obtained in section 2 in a rigorous way. It can also be seen from ( |l4|) that dynamical scaling still 
holds for this model. This is nontrivial, since as pointed out above there are two correlation lengths in the 
problem. 



5 Conclusions 

We have studied a generalised reaction-diffusion process, and the related model of generalized Glauber 
dynamics for a staggered Ising model. We obtained the following results 

i) The dynamical exponent z relating the relaxation time to the equilibrium correlation length is non- 
universal, depending on the ratios of rates or the ratio between the couplings of the Ising model. 

ii) We found that dynamical scaling still holds, i.e. the dispersion relation still depends on a function of 
a single parameter k£ where k is the wave- vector and £ is the equilibrium correlation length. 

It would be interesting to investigate the exact solution of this problem from the point of view of the 
particle dynamics, since one can diagonalize the dynamical operator T by using free fermions. Also it would 
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be interesting to study if there is any change in the dynamical exponent z when one goes away from the free 
fermion condition. One would expect, based on the general argument given in section 2 that the relationship 
of z to the basic rates does not change provided that e'/e 3> pa/pb ^ 1- This means that in the long-time 
limit the system has a very low density of particles and therefore one should expect that the interactions are 



unimportant. Preliminary simulations seem to confirm this result 30 
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